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ANALYTIC APPROXIMATIONS FOR APPLICATION TO SPACECRAFT 
OPTICAL NAVIGATION ON SHORT DATA ARCS 


By Alton P. Mayo and William M. Adams, Jr. 

Langley Research Center 

SUMMARY 

The standard least-squares orbit determination procedures are simplified for 
application to a navigation system where position fixes are taken over a short trajectory 
arc. Basic assumptions used are as follows: (1) for small perturbations in spacecraft 
state, velocity deviations remain constant and position deviations increase linearly with 
time and (2) the sensitivity of the measurements to small spacecraft -position changes 
remains constant over the trajectory arc. Simple analytic expressions are then derived 
for estimating the accuracy of the spacecraft position and velocity as a function of the 
number of fixes and the time span over which the fixes are made, or inversely the num- 
ber of fixes and the time span over which they must be made for a prescribed position 
and velocity accuracy. For a large number of measurements the standard deviation of 
the position error, was shown to vary inversely as the square root of the number of 
measurements, as expected, and the standard deviation of the velocity error was shown 
to vary inversely as the product of the square root of the number of measurements and 
the time span over which the measurements are made. Simplified equations are also 
derived for determining the orbit (state) of the spacecraft, 

INTRODUCTION 

Spacecraft navigation is concerned with both the problem of the single fix, effec- 
tively taking three simultaneous measurements to determine the spacecraft position, and 
with the problem of determining the orbit from numerous measurements taken along the 
trajectory. The use of a single position fix in spacecraft navigation has been the subject 
of several papers (for example, refs. 1 and 2). The development of navigational proce- 
dures for combining multiple fixes or measurements to obtain an estimate of the orbit 
has also been the subject of many papers (for example, refs. 3, 4, and 5). These latter 
procedures are essentially based on modified least- squares techniques and are appli- 
cable to combining measurements made over a trajectory arc of arbitrary length. If the 
trajectory arc is short as in orbit verification calculations or rapid orbit estimation 
techniques, simplifying assumptions may be made in the usual least- squares procedures. 


In the present paper, onboard navigation procedures are developed for measure- 
ments made on a short trajectory arc. The following two simplifying assumptions are 
made in the least- squares procedures: (1) the gravity field is a constant in a small 

region around any point on the nominal trajectory and (2) the mj^nitude and direction of 
the gradient of the measurement remain constant over a short arc of the trajectory. 
Using these assumptions, equations are established for predicting the number of fixes 
as well as the time span, over which they are to be made, for a prescribed position and 
velocity accuracy. Simplified least-squares equations are also derived for determining 
the orbit and estimating its accuracy. 

The procedures presented are useful for the design of onboard orbit determination 
systems using multiple position fixes. The simplified equations presented for determi- 
nation of the orbit from short data arcs are adaptable for either manual or machine 
computation. The equations are illustrated by an example and the results compare 
favorably to those obtained from the complete least- squares solution of reference 6. 

SYMBOLS 


[G] 

h 

I3 

[m] 


matrix relating measurement deviations to perturbations in spacecraft 
position vector 

]l/2 

/a-rr\^ 

magnitude of measurement gradient 


2 

Oz 


0yr 

dx) ""\9y 


unit matrix of rank 3 


matrix relating measurement deviations to perturbations in the space- 
craft state vector 


N 

[p] 


T 

[w] 


number of time points when measurements are made, that is, number of 
fixes 

null matrix of 3 X 3 dimensions 

local position information matrix; inverse of position covariance matrix 
time from epoch 
length of data arc 

measurement weighting matrix, a diagonal of 


2 



! 


x,y,z perturbations in spacecraft geocentric equatorial coordinates 

y deviation in optical measurement 

Y vector of coefficients in least- squares normal equations 


€ 

random error in optical measurement 


O’ 

standard deviation 



standard deviation of optical measurement 


Or 

square root of trace of position covariance matrix, 


CTv 

square root of trace of velocity covariance matrix, 


[^] 

transition matrix 


« 

column matrix of the state vector 


{y> 

column matrix of observations 


Subscripts: 



i 

summation index 


o 

at epoch 


P 

position 


t 

time from epoch 


to 

at time of the epoch 


x,y,z 

derivative with respect to x,y,z, respectively 



derivative with respect to x,y,z, respectively 
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Matrix notations: 


{ > 


column matrix 


n 

fw] 

(info] 

[]T 


square or rectangular matrix 


covariance matrix 


information matrix (inverse of covariance matrix) 


transpose of matrix {] ] 


Optical measurement symbols: 

subtended angle of Moon 

observation of an angle from a star to Earth horizon 




* ^ 

+ 1 ) 

© 

© 

0 

o> 


observation of an angle from a star to a Moon landmark 


observation of an angle from a star to Moon horizon 


observation of an angle from a star to an Earth landmark 


observation of Earth azimuth 


observation of Earth elevation 


angular diameter of Earth 


^ -He „ observation to star 1, star 2, and star 3, respectively 
If 2> 3 

Dots over symbols denote derivatives with respect to time. 


WEIGHTED LEAST- SQUARES PROCEDURES 


The weighted least-squares procedures for spacecraft navigation are based on the 
assumption that small trajectory perturbations are transitioned linearly along the tra- 
jectory. This concept may be expressed mathematically as 
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( 1 ) 


{^> = W <=Ho> 

where the column matrices and are the trajectory perturbations at 

time t and t^, respectively. 

The elements of the so-called transition matrix |0] are partial derivatives 
which relate perturbations in spacecraft state at time t due to perturbations in state 
at time t^. The perturbations in the observations made on the trajectory are also 
assumed to be linearly related to the trajectory perturbations; this leads to the 
e^spression 


(y) =[m]{Xj) + {€} (2) 

The column matrix {y} is the difference between the measurement and the 
values calculated for the estimated trajectory. The state vector {Xj.^ represents the 
deviation of the actual trajectory from the estimated trajectory. 

Combining equations (1) and (2) gives equations of condition, which in matrix form 
are 


{y> = MM{xto> + (0 (3) 

It is the usual procedure to multiply each observation equation by the reciprocal 
of the measuring accuracy; this gives large multiplying or weightii^ numbers to those 
equations associated with high measuring accuracy and small weighting numbers to those 
equations with low measuring accuracy. Thus, the observation equations become 

These are a set of linear equations expressing the deviations in the observations 

at time t due to perturbations in the spacecraft state at time t^ and the random error 

introduced when the observations were made. Given a set of observational data and 

estimated values of the observations, the weighted least-squares technique may be used 

to determine the trajectory perturbation This involves the fitting of equation (4) 

minus the error term to the data. Multiplying the modified equation (4) by 
T 

and solving the resulting equation by matrix inversion yield the 
weighted least- squares solution as 
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For normally distributed, uncorrelated measurement errors, the matrix 


[MT[„]T[w][m]Mj' (6) 

is the covariance matrix of the estimated and expresses the accuracy of this 

estimate from the knowledge of the measurement accuracy as expressed by the 
matrix [w^ and by using the linear relation expressed by [m] and 

In subsequent sections of this paper, it is shown that the matrix [m^ can be 
assumed to be constant over a short trajectory arc and that the matrix [0^ can be 
simply approximated by assumir^ linear perturbative motion. This, in addition to sim- 
plifying the equations, yields insight into the behavior of onboard navigation procedures. 

APPROXIMATION TO THE TRANSITION MATRIX 


For the weak gravitational field in most of translunar and interplanetary space, the 
approximation can be made that x^ = x^ = x^ = yy = y^ = Zjj " ^y “ ^ ^ small 

region about any point on the nominal trajectory. Then the equations, for the perturbative 
motion, reduce to 


^t ■ ^o 

^t = ^o / 


(7) 


From which it follows that 


Xj = Xq + 

n = ^o yo* 

^t = S + ^o* 

and the perturbative motions are linear. In this case, the transition matrix [0^ 
then be expressed as: 



I 


1 0 0 t 0 0 


0 1 0 0 t 0 


m- 


0 0 0 


0 

1 


0 t 
0 0 


(8) 


0 

10 


0 

0 


0 

0 


0 

0 


1 

0 


0 

1 


The transition matrices obtained from integration of the complete equations of 
motion, that is, with and yy and so forth not equal to zero, are shown in table I for 
a translunar trajectory. Selected elements of the transition matrices obtained from inte- 
grating the complete equations of motion are also plotted in figure 1. By examining equa- 
tion (8), figure 1, and table 1, it can be seen that the approximation of linear perturbative 
motion is good except in the near-earth region (flight time = 0). When comparir^ the 
various curves of figure 1, note that the linear position variation, for 10, 20, and 52 hours 
of flight time, is masked by the highly sensitive vertical scale. In any near -body region 
a more complete transition matrix could be employed such as the ejq)ressions for the 
transition matrix given in reference 7. The diagonal elements of the transition matrix 
for values of t up to 8 days on a Mars approach trajectory are shown in figure 2. 

APPROXIMATION TO THE GRADIENTS OF THE MEASUREMENTS 

The relation between the local measurements and the spacecraft local position is 
given by 




X 


^1 


II 

y 


^2 ^ 

^3 
^ J 


z 

•k. -j 


^3 
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I I 


L 


II 


where 


[G] = 


9yj 

9yi 

9yi 

dx 


9Z 




dx 

9y 

9Z 


ayg 


dx 

9y 

9z 


For optical measurements, the 3x3 matrix [gJ expresses the gradients of the 
measurements with spacecraft position. For short data arcs, say, 4 hours of a trans- 
lunar orbit, the direction and magnitude of the gradient vector should not change con- 
siderably, except in the regions where the spacecraft is near the body on which the meas- 
urements are being made. A plot of the magnitude of the gradient of the measurements 
over a 4 -hour interval for flight in various portions of translunar space is shown in fig- 
ure 3 which indicates that the magnitudes are fairly constant except for a measurement 
made near the earth at 10 hours and measurements made near the moon at 52 hours. 
Thus, if one avoids the near -body regions, the gradients should remain approximately 
constant over a 4 -hour translunar flight time or, in order to maintain approximately con- 
stant gradient over the arc, the length of the data arc has to be shortened in the near- 
body region. 

The results for the Mars mission are shown in figure 4, The plot shows that the 
magnitudes of the gradients are fairly constant over 1-day segments in the near -body 
regions cf the Mars trajectory. The magnitude of gradients would be expected to be 
more nearly constant in the far -body region. 

DETERMINATION OF THE STATISTICS OF THE ORBIT 

Accuracy estimates of the spacecraft state are a necessary part of navigational 
computations. These estimates are usually in the form of a covariance matrix whose 
diagonal elements are the square of the 1 -sigma uncertainty of the spacecraft state 
vector. A simplified expression for the covariance matrix is derived by usir^ the pre- 
viously mentioned assumptions of a short data arc, linear propagation of small trajectory 
perturbations, and constant sensitivities of the optical measurements with respect to 
small deviations in spacecraft position. 

The equation relating the spacecraft position to the measurements is given by equa- 
tion (9) and in referring to the reasoning in equations (4) to (6) the associated covariance 
matrix of position is 
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( 10 ) 



-1 


where W is the measurement weighting matrix. The local state information matrix 
due to the local measurement is thus given by 


{inf^ 


t " 


g'^WG ' Og 
1 


* 

O 

1 

1 J 

1 

_ *^3 ' ®3 _ 


V 

jz : »3_ 


( 11 ) 


where G^WG is the upper 3x3crfa6x6 matrix and is designated as P. The infor- 
mation matrix of the state at time due to the local fix (three simultaneous observa- 
tions) is given by 

[Mo],o= 


which, by using equation (8), becomes 


[Inf^to = 


lie 


L«3 


tio 


t^I. 


( 12 ) 


and the information matrix at time t^ due to all fixes over the data interval is given by 
the summation process 


1 

N 

Nig : 

j 

1 

t 

i=l 


N 

7 — ' 

V ' 

z *1%: 

Z 

i"i , 

i=l 




( 13 ) 


where, as in equation (13), P is assumed constant for all data points and N, 
N 


X 



are scalar coefficients of 3 X 3 identity matrices. 


This expression may 


i=l 

written in closed form as 



i=l 

be 


and 


NI, 


[inf^ 


to 


NT T 


NT , 
“^3 


(14) 


N(2N - 1)T^ j 
6(N - 1) 3 


where the assumptions have been made that the first data point is obtained at t = t^ = 0.0 
and the spacing between data points is constant. 

The covariance matrix associated with the estimate of x^^ is given by 


-1 


[covx^ =[lnf^^o = 


2(2N - 1) j 
N(N + 1) 3 


-6(N - 1) j 
N(N + 1)T 3 


-6(N X 1) j 
N(N + 1)T 3 


.-1 


(15) 


N(N + 1)T^ 


.-1 


A comparison of the results of this equation and the results of reference 6 for a 
navigation system utilizing 4 hours of data beginning at 20 hours of flight time is shown 
in table II. It is seen that the results from the approximate solution compare favorably 
to those obtained from the complete least-squares solution. 

Examination of equation (15) shows that the relation between the position and veloc- 
ity uncertainty at any time t is given by 

^ = T ^ for large N (16) 

The comparison of this equation with the complete least -squares results of reference 6 
is given in figure 5, In this figure the data span is varied from 0 to 4 hours in various 
translunar regions. As can be seen, the agreement is good from 10 to 40 hours of flight 
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time. It should be stated here that the ratio of the position to the velocity uncertainty is 
not a function of the type of measurements being made. 

This same position -to -velocity uncertainty relation is shown in figure 6 for a 
4-hour data span at various distances from earth for numerous navigation systems. 

VARIATION OF POSITION AND VELOCITY ACCURACY WITH 

NUMBER OF FIXES 


From equation (15), it is seen that 


2 . 2(2N - 1) 
r N(N + 1) 


• Trace [^P]] 


-1 


Prom which it follows that 


cr 


r 


2(2N - 1) 

N(N + 1) r,o 


for large N 

s/n 


Also from equation (15) 


(17) 


or 



12(N - 1) . Trace [p]“l 
N(N + 1)T^ 


= 2 



N - 1 
N(N + 1) 


^r.o, 2\[3 
T 


for large 


N 


(18) 


Plots showing a comparison of results from equations (17) and (18) with the complete 
least-squares results of reference 6 are given in figures 7 to 10. The values used for 
CFj. ^ in figures 7, 8, and 9 are shown in figure 10. This figure shows a comparison of 
the values of a single position fix imcertainty <r obtained from equation (17) and 
those obtained from diagonal elements of the equation 

[covx„]= [oqviog] 


(19) 



DETERMINATION OF NUMBER OF FIXES AND TIME SPAN TO OBTAIN 
A PRESCRIBED POSITION AND VELOCITY ACCURACY 


The niimber of measurements to obtain a prescribed by using a navigation 
system with accuracy a for a single fix, is given from equation (17) as 



For large values of N, equation (17) reduces to 


Thus 




for large 


N 


The time span over which the measurements must be made to obtain a specified 
velocity uncertainty along with this 0 ^, is given from equation (16) as 

O'- l~ 

— as — ^3 for large N (21) 

Uy (Ty 



USE OF LESS THAN THREE OBSERVATIONS PER TIME POINT 


From equation (13), the information matrix is approximated by the equation 



This equation does not require that three measurements (one fix) be made per time 
point. The G matrix could have the second or second and third row zeros. This 
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II 


would correspond to making one or two measurements per time point. However, in this 
case, the information matrix would not have an inverse since the P matrix of the above 
product has no inverse. Thus the simplified expressions relating Oj, and Oy in terms 
of (Tj. p, N, and T developed in this paper are only applicable when three or more 
measurements are made per time point. Nevertheless, if a priori information is avail- 
able, it is possible to obtain an improved estimate of the orbit with less than three obser- 
vations per time point while retaining the simplifying assumptions about the gradients of 
the measurements and the time variation of the transition matrices. A comparison of 
some results from reference 6 shows that in the unrestricted orbit determination pro- 
cess, the three observations per time point cases give considerably better results than 
the two observations per time point cases; however, the orbit was determined in both 
cases. These comparisons are shown in figure 11. 

DETERMINATION OF THE SPACECRAFT STATE USING 
THE APPROXIMATE TRANSITION MATRIX 


The weighted least-squares normal equations (eq. (5)) for the spacecraft state are 
of the form 

1-1 


{^o} = ^o]' 

where the column is given by 


(y}= [ 0 ]T[m]T[w](y}= ^ ^ 



(23) 


(24) 
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From equations (8) and (9), is assumed constant over the data interval and 


9Yi 

— - is assumed given by 

9*0 


9X 


9Yi 9Yi 


(25) 


thus, the expression for can be written from equations (24), (25), and (9) as 



N * 

i=l^^i 
N ^ 

i=l ^i 
N - 

Y ! 3 i 

^ 0.2 

i=l *^i 

Y yii^ 

Z/ rr 2 




i=l ‘^i* 


N 


til V 

Y ^3i‘i 

.A 0.2 

1=1 1 


(26) 


Combining this equation with equation (15) yields the solution for spacecraft state as 



2(2N - 1) j 

6(N-1) J 


P"^ 1 Oo 


N(N + 1) ^ 

N(N + 1)T ^ 


' , 3 

{--}= 




» 


6(N - 1) j 
N(N + 1)T 3 

12(N - 1) I 


0« ' 


N(N + 1)T^ 


1 


{v> 


(27) 


A comparison of approximate navigational results obtained from equation (27) with those 
from the complete least-squares solution of reference 6 is shown in table in . The data 
in table III represent 40 fixes made from 20 to 24 hours of flight time. 
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CONCLUDING REMARKS 


Optical navigational techniques applicable to position fixii^ on short trajectory arcs 
were developed. The gradients of the optical measurements over a short arc were 
assumed to be constant in magnitude and direction. Also, It was shown that the gradient 
of the gravity field in a limited region about the nominal trajectory could be assumed to 
be zero. By combining these two assumptions, simplified least-squares normal equations 
for the determination of the orbit and its accuracy were developed. These equations show 
that the standard deviation of the position error determined from a large number of mea- 
surements varies inversely as the square root of the number of measurements, as 
expected, ard the standard deviation of the velocity error varies inversely as the product 
of the square root of the number of measurements and the time span over which the 
measurements are made, From these relations, the time span of the measurements and 
the number of measurements which must be made for a prescribed spacecraft position 
and velocity accuracy was readily established. The approximate least-squares equations 
for determining the orbit were applied to an example and the results compared favorably 
to those obtained from the complete least -squares solutions. 

Lar^ley Research Center, 

National Aeronautics and Space Administration, 

Lar^ley Station, Hampton, Va., June 11, 1968, 

125-17-05-09-23. 
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TABLE I.- TRANSITION MATRICES AT 20 AND 52 HOURS OF FLIGHT TIME 


IN A TRANSLUNAR TRAJECTORY 


Transition matrix | 

from 20 to 23.9 hours of flight time: 




0.9945 

0.29 X 10-3 

0.19 X 10-3 

14015 

1.1 

0.58 


.98 X 10-3 

1.0078 

.68 X 10-2 

1.0 

14070 

31 


.0 

.63 X 10-2 

.9980 

.52 

31 

14032 


-.74 X 10"® 

.37 X 10-*? 

.20 X 10-"? 

.9950 

.18 X 10-3 

.93 X 10- 


.37 X 10-’^ 

.98 X 10-8 

.96 X 10-6 

.18 X 10-3 

1.0064 

.63X 10- 


.19 X 10-'? 

.95 XlO-8 

-.22 X 10-8 

.92 X 10-4 

.63 X 10-2 

.9885 

Transition matrix | 

from 52 to 

55.9 hours of flight time: 




0.9997 

0.15 X 10"3 

0.24 X 10-3 

14040 

7.1 

3.2 


.0 

.9999 

.39 X 10-2 

6.2 

14045 

9.3 


.97 X 10-3 

.20 X 10-2 

.9990 

3.6 

9.8 

14033 


.75 X 10-8 

.22 X 10-8 

.10 X 10-8 

1.0002 

.17 X 10'2 

.83 X 10- 


.22 X 10-8 

.19 X 10-8 

.29 X 10"® 

.17 X 10-2 

1.0013 

.21 X 10- 


.10 X 10-8 

.29 X 10-8 

-.20 X 10-8 

.83 X 10-3 

.21 X 10-2 

.9984 


TABLE n.- COMPARISON OF COMPUTED COVARIANCE MATRIX WITH 
ANALYTIC APPROXIMATION TO COVARIANCE MATRIX 

[Sitnilar to Right time of 20 hours of ref. 6; 100 fixes at 1 per 
0.04 hour; star to earth horizon, two stars to moon horizo:^ 


Computed least-squares values; 


COV XQ = 


5.2876 

5.8115 

4.0730 

-0.5675 X 10-3 

-0.6031 X 10-3 

-0.4211 X 10-3 

5.8115 

22.1179 

12,2335 

-.6077 X 10-3 

-2.2777 X 10"3 

-1.2597 X 10-3 

4.0730 

12.2335 

10.0330 

-.4232 X 10-3 

-1.2573 X 10-3 

-1.0323 X 10-3 

-.5675 X 10-3 

-.6031 X 10-3 

-.4211 X 10-3 

.8088 X lO-"^ 

.8379 X 10-’^ 

.5792 X 10-'? 

-.6077 X 10-3 

-2.2777 X 10-3 

-1.2596 X 10-3 

.8379 X lO-*^ 

.3063 X 10-3 

1.6888 X 10-*? 

-.4232 X 10-3 

-1.2573 X 10-3 

-1.0323 X 10-3 

.5792 X lO""^ 

1.6888 X 10-’? 

1.3913 X lO-"? 


Analytic approximate values: 



5.3431 

6.0529 

4.2825 

-0.5565 X 10-3 

-0.6305 X 10-3 

-0.4461 X 10“3 

6.0529 

24.0818 

13.3539 

-.6305 X 10-3 

-2.5085 X 10-3 

-1.3910 X 10-3 

4.2825 

13.3539 

10.8570 

-.4461 X 10-3 

-1.3910 X 10“3 

-1.1309 X 10-3 

-.5565 X 10-3 

-.6305 X 10-3 

-.4461 X 10-3 

.7703 X lO-*? 

.8757 X lO-”? 

,6195 X 10- ■? 

-.6305 X 10-3 

-2.5085 X 10-3 

-1.3910 X 10-3 

.8757 X 10-"? 

.3484 X 10-3 

1.9319 X lO-"? 

-.4461 X 10-3 

-1.3910 X 10-3 

-1.1309 X 10-3 

.6195 X 10-"? 

1.9319 X lO-*? 

1.5707 X lO-*? 


TABLE in.- EXAMPLE OF APPROXIMATE AND COMPLETE 
LEAST-SQUARES SOLUTION FOR SPACECRAFT STATE 
ON A TRANSLUNAR TRAJECTORY 


Parameter 

True deviation 
of state from 
nominal value 

Least -squares 
solution of 
reference 6 

Approximate solution 
of this paper 

X, km 

9.656 

11.4 

8.5 

y, km 

9.656 

9.6 

5.5 

z, km 

9.656 

15.4 

10.8 

X, km/ sec 

0.0061 

0.0062 

0.0066 

y, km/sec 

0.0061 

0.0062 

0.0072 

z, km/sec 

0.0061 

0.0055 

0.0067 


0 .8 1.6 2,h 3*2 ‘^-.C 

Time since epoch^ hr 

(al Position, 


Figure 1.- Trajectory displacements from nominal due to un 
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1.00 

. 9 ? 




Figure 2.- Variation of elements of state transition matrix over an 8-day span near Mars. Range to Mars from 3 270 000 tern to 289 000 km. 
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Description of observations 
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Figure 6.- Comparison of simplified expressions for Oiyby with frans lunar data from reference 6, N = 40; T = 14 040 sec* 




o> 


Description of observation 





(a) Constant observing time (2 hrL 


(b) Constant number of fixes (401 


Figure 7*- Comparison of simplified expression for navigation accuracy with numerical results of reference 6* 
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